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ABSTRACT 

While the nature of the heating mechanisms in the corona remains elusive 
their associated cooling is also a poorly known but a far less observationally re¬ 
strictive subject. In this work, we analyse coordinated observations spanning 
chromospheric, transition region (TR) and coronal temperatures at very high 
resolution which reveal essential characteristics of thermally unstable plasmas. 
Coronal rain is found to be a highly multi-thermal phenomenon with a high 
degree of co-spatiality in the multi-wavelength emission. EUV darkening and 
quasi-periodic intensity variations are found to be strongly correlated to coronal 
rain and especially ‘showers’. Progressive cooling of coronal rain is observed, 
leading to a height dependence of the emission. Furthermore, a fast-slow two- 
step catastrophic cooling progression is found, which may reflect the transition 
to optically thick plasma states. The intermittent and clumpy appearance of 
coronal rain at coronal heights becomes more continuous and persistent at chro¬ 
mospheric heights just before impact, in agreement with previous observations 
above sunspots. This change of character is mainly due to a funnel effect from the 
observed expansion of the magnetic held at low heights. Strong density inhomo¬ 
geneities on spatial scales of 0.2" — 0.5" are found, in which TR to chromospheric 
temperature transition occurs at the lowest detectable scales. The shape of the 
distribution of coronal rain widths is found to be independent of temperature 
with peaks close to the resolution limit of each telescope, ranging from 0.2" to 
0.8". However we hnd a sharp increase of clump numbers at the coolest wave¬ 
lengths and especially at higher resolution, suggesting that the bulk of the rain 
distribution remains undetected. Rain clumps appear organised in strands. Such 
structure is not limited to chromospheric temperatures but extends at least to 
TR temperatures as well, suggesting an important role of thermal instability in 
the shaping of fundamental loop substructure. At the smallest detected scales 
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we further find structure reminiscent of the MHD thermal mode. Rain core den¬ 
sities are estimated to vary between 2 x 10^° cm“^ and 2.5 x 10^^ cm“^ leading 
to signihcant downward mass fluxes per loop of 1 — 5 x 10® g s“^, suggesting a 
major role in the chromosphere-corona mass cycle. 

Subject headings: magnetohydrodynamics (MHD) — Sun: activity — Sun: corona 
— Sun: hlaments, prominences 


1. Introduction 


The solar corona, the outermost layer of the solar atmosphere, has an average tem¬ 
perature above a million degrees, a fact that has puzzled astrophysicists for decades. The 
processes responsible for this heating are still unknown, mostly due to the very short time 
and spatial scales on which they operate. However, not everything in the corona is hot. 
Structures at chromospheric temperatures are clearly present in the solar corona. Evidence 
for this has been found for more tha n a century, since the hrst detection of prominences 
in the solar atmosphere (jSecchilll875l) . In the last 40 years, especially in the last few years 
with the advent of high resolution instruments, smaller chromospheric structures with a 
more pervasi ve character hav e been shown to occur intermittently in the corona above ac¬ 


tive regions flKawaguchi 
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Ahn et al 


20141) . This material is known as coronal rain, partially 


ionised, dense clumpy structures catastrophically cooling in the corona and accreting to¬ 
wards the solar surface. Its occurrence is considered as the direct observational consequence 
of the thermal instability mechanism, in which radiative losses locally overcome the heating 
processes. 


The presence of prominences in the corona and especially the recent detection of coro¬ 
nal rain with a more pervasive character entails fundamental aspects of the interaction of 
plasmas and helds in low f3 environment in which thermal instability plays a central role. 
This mechanism is intrinsically linked to the characteristic cooling aspects of the heating 
mechanisms and provides far less strict detection possibilities due to the longer timescales 
and an enormous gain in spatial resolution when observing in wavelengths associated with 
cool emission. The study of coronal heating through associated cooling therefore provides 
unique advantages. 


A big unknown in the physics of the solar corona is what percentage of coronal plasma 
is thermally unstable. This is also of signihcant importance for stellar coronae, since high 
speed accretion of chromospheric plasma can produce characteristic prohles in their UV and 
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optical ran ge spectra, both from the falling plasma and from the impact into the lower at¬ 
mosphere f Avres fc Francell2010l: iReale et ahll2013l. 120141) . Based on numerical simulations, 
such conditions me expected if the heating is mainly conce ntrated towards the footpoints 
flAntiochos et alJll999l: iKarpen et al.ll200ll: iLuna et al.ll2012l) . a condition often reported in 
active regions, and are to a lesser ext ent related to the tempor al character of the heating, 
i.e. whether it is impulsive or steady (iMiiller et al.ll2003L 120041) . Forward modelling of nu¬ 
merical results and comparison with EUV light curves obtained in observations have helped 
answering this question, but a debate exists on whether the heating requirements setting 
loops in a thermal non-equilibrium state fully match the g e neral observed characteristics of 


activ e region loops fiMok et ahl 120081: iKlimchuk et al.l 120101 iPeter et al.l 120121: iLionello et ah 
20131: iKlimchukl 120141) . Coronal rain observations can help answer this question by esti¬ 


mating the coron al volume and m ass flux involved in these events. A Erst attempt at this 
was performed by ISchriiverl (120011) using channels from the Transition and Coronal Explorer 
(TRACE), which lead to an occurrence interval of several days for one active region loop to 
undergo catastrophic coolin g. Based on Ha data from the CRI S P ins trument at the Swedish 
1-m Solar Telescope (SST). lAntolin fc Rouppe van der VoortI (120121. hereafter Paper 1) re¬ 
duced the occurrence frequency to half a day or less, and estimated the coronal rain mass 
flux to be about one third of that uploaded into the corona from spicules, thus stressing an 
important role in the chromosphere-corona mass cycle. 

Paper 1 predicted that most of the rain material may pass undetected in coarse reso¬ 
lution instruments due to average clump0 widths of 300 km (and lengths of 700 km or so). 
Condensations were reported to be clustered in space into groups of correlated events, which 
indicates that the properties of coronal heating can be correlated over a signihcant spatial 
scale and that individual strands (the theoretical internal components of coronal loops) are 
not necessarily uncorrelated in their evolution. However, while most condensations were 
found to be correlated in time, only few events, termed ‘showers’ in Paper 1, were found 
to occur close enough, which would produce large enough EUV absorption detectable with 
the coarser resolution of coronal instruments. Still, the limit cycles that thermally unstable 
loops experience, predicted from numerical simulations, should be visible observationally as 
periodic variations of the EUV intensity (which would then be correlated to the appearance 
of coronal rain). Such behaviour has recently been rep orted and seems to £t well into the 
thermal non-equilibrium scenario (IFroment et al.l 120141) . 


Since it is chromospheric (cool) plasma occurring in the corona, coronal rain constitutes 
the highest resolution window into the substructure of loops, and provides a tool for reveal- 
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ing the global magnetic field topology. First evidence from this was provided in Paper 1, 
where tracking of rain over a decaying active region allowed us to calculate the fall angle 
of the clumps, and therefore to reveal the large scale morphology of the magnetic field and 
distinguish several loop families by tracing the clumps along the legs. The coronal rain con¬ 
densations were reported to be organised locally as strand-like structure with average widths 
around 300 km and lengths around 700 km, with nonetheless broad distributions down to 
100 km and 200 km, and up to 700 km and 2400 km respectively for widths and lengths. 
The multi-strand su bstructure of toops w as also reported in on- disc observations of c oronal 
rain with CRISP in lAntolin et al.l (120121) and more recently in I Scullion et al.l (120141) . The 
latter work extends such findings to the coronal rain produced in post-flare loops and to the 
apexes of loops. In that work the distribution of widths peaks at 100 km, clearly indicating 
that the bulk of the distribution is unresolved due to the lack of resolution. Interestingly, 
the lengths of some of the observed condensations were reported to span up to 26000 km, 
further underlining the global coronal magnetic field tracing possibilities offered by coronal 


ram. 


The strand-like structure revealed by coronal rain along loops appears similar to the 
thread-l ike structure observe d in prominences of certain types (e.g., active region promi¬ 


nences, lOkamoto et al.N2007l) . A relevant question is therefore whether such structure is 
strongly dependent on thermal instability, or whether such structure also exists at coronal 
temperatures und er equilibrium conditions, in accordance with the debated multi-strand 
scenario in loops (Klimchuk 2006 : Brooks et ah 2012 . 2013 : Peter et ah 2013). The loss of 
pressure produced by thermal instability produces local accumulation of material leading to 
condensations. Given a thermally unstable loop, the continuous occurrence in the same spot 
of thermal instability, combined with the produced strong flows (or other possible siphon 
flows) in a low (3 environment such as the cor ona can easily lea d to strand-like structure, as 
shown in multidimensional MHD simulations (IFang et al.ll2013j) . In this paper we show that 
even in this scenario it is reasonable to expect that such strand-like structure is not locally 
confined to the thermally unstable regions (and thus only confined to the low temperature 
plasma in the loop), but extends along the entire loop and therefore to the hot coronal 
plasma as well. Therefore, even if a loop does not tend to be originally multi-stranded, it is 
possible that such organisation is attained due to thermal instability and maintained there- 
afte r. Such scenario, for which we provide support in this paper (also theoretically supported 
by ILow et al.ll2012bl) . has not been considered before and may significantly impact on the 
evolution of the coronal loop. 

It can therefore be expected that a multi-thermal structure accompanies a multi-strand 
organ isation of the plasma in thermally unstable loops. This is shown in IScullion et ah 


(120141) , where loops observed in the coronal filters of AIA are co-spatial with the Ha conden- 
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sations. The lack of resolution of AIA unfortunately does not allow us to see the details and 
spatial scales at which such multi-thermality can exist. In this work, with the help of the 
high resolution of the Interface Region Imaging Spectrograph (IRIS), we extend this result. 

We present two datasets that combine multi-wavelength instruments spanning chromo¬ 
spheric, transition region and coronal temperatures. SST/CRISP , Hinode/SOT, IRIS/SR 
and SDO/AIA are used to reveal the multi-strand and multi-thermal aspects of coronal rain 
at unprecedented detail. The paper is organised as follows. In Section [2] the observations 
are presented. In Sections [3] and 0] the temperature and morphology of the observed conden¬ 
sations are presented, respectively. Results are discussed in Section |5] and conclusions are 
given in Section |6l 


2. Observations 


2.1. Data processing 


In this work we analyse two different datasets corresponding to solar limb observations 
above active regions. The hrs t dataset combines data from the Swedish 1-m Solar Telescope 
fS'S'T: IScharrner et al.ll2003al) with the CRISP spectropolarimeter and the AIA instrument 
( Lemen et ah 20 111) of the Solar Dynamics Observatory (SDO), dates from 26 June 2010 from 
10:03 UT to 11:40 UT, and was centred on AR 11084 at b,?/] = [-875, -319 ]. The second 
dataset combines the SOT telescope on board of Hinode (iTsuneta et al.ll2008l) . the Slit-Jaw 
Imaging (SJI) instrument of IRIS I De Pontieu et ah 2014 ) and SDO/AIA, dates from 29 
November 2013 from 22:30 UT to 23:30 UT and focused on AR 11903. IRIS/SR was centred 
at [x, y] = [944, —264] fully containing the Hinode/SOT FOV, centred at [x, y] = [959, —220]. 
Context images for these coordinated observations are presented in Figs. [T] and [2l where 
different colours denote different wavelengths and the different held-of-views (FOV) of each 
instrument can be compared. 


We mainly focused on the SDO/AIA hlters 304, 171 and 193, corresponding to a range 
of transition region to coronal temperatures usually spanned by thermally unstable plasmas. 
The upper photospheric hlters 1600 and 1700 were also used for co-alignment purposes (see 
Section [2.2p . Level 1.5 data from the SDO/AIA was used, obtained through the normal 
calibration routines in SolarSoft. Full disc images are provided from AIA with a cadence of 
12 s, an image scale pixel of 076 pixeU^ and a spatial resolution of 173 - 177. 

IRIS combines a Slit-Jaw imager (SJI) and a spectrograph (SG). The IRIS data used 
here comprises only the SJI instrument since the IRIS slit was not located on the region 
of interest. Therefore, in this work only the coronal rain imaging capabilities of IRIS are 
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shown. The dataset corresponds to level 2 data from a sit-and-stare observation combining 
the SJI 1400, 1330 and 2796 hlters, dominated respectively by the Si IV transition region 
lines 1393.78 A and 1402.77 A (formed around logT = 4.8), the C II chromospheric lines 
1334.53 A and 1335.71 A (formed around logT = 4.3), and the Mg II k chromospheric line 
2796.35 A (formed around logT = 4). IRIS/SJl provides a FOV of 175" x 175", an image 
scale pixel of 07166 pixel“^ and a spatial resolution of 0733 (for the FUV) - 074 (for the 
NUV). The cadence for this observation is 36.5 s in each hlter. 


The CRisp Imaging SpectroPolarimeter fCRISP nScharmer et al.ll200811 at the SST sam¬ 
pled the Ha spectral line at 41 line positions with 0.085 A steps (3.9 km s“^), in the spectral 
window [—1.7,1.7] A (~ ±78 km s“^) with a cadence of 9 s. The FOV of the SST for 
the pres ent observation is 53" x 57". The images benehted from the SST adaptive optics 
system fjScharmer et ahl l2003bl) and the image restoration tec hnique Multi-Object Multi- 
Frame Blind Deconvolution iMOMFBD dvan Noort et al.ll2005h . Although the observations 
suffered from seeing effects, most of the images are close to the theoretical diffraction limit 
for the SST at the wavelength of Ha: A/D ~ 0714. We employed the same reduction pro- 
cedure as in Paper 1 which used early versions of parts of the data pipeline CRISPRED 
(fdc la Cruz Rodriguez et~^ 2015 ). 


Hinode/SOT recorded hltergrams in the Ca II H band (formed around logT = 4) at 
a cadence of 4.8 s. The image scale pixel is 07109 pixeR^ and the spatial resolution is 
A/D ~ 072. The SOT FOV for this observation is 55.8" x 55.8". Processing of data was 
carried out through the normal calibration routines in SolarSoft. 


2.2. Co-alignment procedures 

The co-alignment of S'S'T/CRISP and SDO/AIA data for the 26 June 2010 hrst con¬ 
sisted in determining the solax-x/y coordinates and locate the solar north direction for each 
CRISP frame. We then eliminated the SST drift by hxing each frame on the sunspot of 
AR 11084, at the centre of the SST FOV, to reduce the frame-to-frame motion of the loop’s 
footpoints (and therefore of the loops themselves). The location of the solar limb for each 
frame was also determined, which allowed to apply a gradient hlter reducing the intensity 
of on-disc features for better visualisation of the off-limb structures. SDO data was then 
cropped and interpolated to the exact SST FOV. Small misalignments still remain, which are 
then corrected through cross-correlation using photospheric bright points (with long enough 
lifetime) and the solar limb location (using the AIA hlters 1600 and 1700, in combination 
with the far wing positions of Ha, which mostly show the photosphere). 
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For the SOT-SJI-AIA co-alignment, we first removed the slow drift of SOT by co- 
alignment of the images through tracking the movement of the limb. The cumulative offsets 
from this drift were applied to the time series, resulting in a rigid co-alignment on the order 
of one pixel or less. The IRIS/SJl images do not have this slow drift, but occasional wobble 
shifted the images a few pixels (this is normal during eclipse season). Such shifts were 
measured by applying cross-correlation techniques and applied to the time series. Then, to 
co-align the images from different instruments the following procedure was used. First, the 
AIA images were co-aligned and cropped to the SJI held of view by identifying common 
bright points in the AIA 1600 hlter and the SJI 1400 hlter. Afterwards, both AIA and 
SJI images were interpolated to SOT’s image pixel scale and were co-aligned using common 
features on the SJI 1400 hlter and SOT’s Ca II H hlter. 


2.3. Methods 


Due to the complex interplay of forces in loops coronal rain dynamics can be extremely 
complex, especially when observing at sub-arcsecond resolution. Furthermore, their thermo¬ 
dynamic sta te is also expected to change, as we show in this paper, and as has been also 
indicated in lUarra et al.l (120141) . A consequence of this is that it is very difficult to follow 


a single rain clump (at high resolution) over long distances. The clumps will usually split, 
merge, elongate along their paths. For these reasons, especially for the CRISP data where 
myriads of small clumps are observed and for which the line-of-sight (LOS) superposition 
is signihcant, individual tracking is very difficult and time consuming. For automatically 
detecting rain in the CRISP data we have therefore opted for perpendicular cuts across the 
loops at positions where rain is observed, rather than longitudinal paths along individual 
clumps. Detection is done through semi-automated procedures verifying specihc conditions 
for the clumps. These consist especially of contrast thresholds (a rain pixel needs to be 
signihcantly brighter or darker with respect to the background, depending on whether the 
detection is realised off-limb or on-disc, respectively; all frames with bad seeing are then 
removed), intensity variation thresholds (a rain pixel must present an average change in 
intensity above a specihc threshold within a reasonably small time interval), and size re¬ 
strictions (a clump must span at least 3 rain pixels in the direction perpendicular to the 
direction of propagation). Once the rain pixels are identihed the clump sizes are calcu¬ 
lated with Gaussian hts along and across the direction of propagation. Therefore, one size 
value for a clump is the average over several tens or hundreds of measurements. The results 
are also tested for specihc conditions, such as size thresholds, small ht errors, consistent 
widths (a clump’s width must not diher signihcantly over a short distance along the axis 
of propagation) and so on. For rain detection with CRISP we have further opted for us- 





ing Doppler intensities rather than the clnmp intensities alone, since detection is improved 
in this way. The Doppler intensity for a rain pixel at position x, y and time t is dehned as 

Aq —(A —Ao), t), where Aq is the theoretical wavelength 

valne at rest. 


For the rain observed with SOT, S JI and AIA we have opted for the ‘traditional’ method 
nsed in previous studies. First, individual clumps’ paths are traced manually. For each time 
step at which rain appears along the path and at each position along the clumps’ paths 
semi-automated routines measure the width and length following the same criteria specihed 
above. The average and standard deviation of these measurements for a given path provide 
the clump’s size (width and length) and corresponding errors. 


For all the above measurements, regardless of the instrument, clump intensities are 
calculated by subtracting from a given rain pixel intensity the average background intensity 
at the same pixel position. The average is calculated over a time interval close to that 
of detection and for which no rain is detected. The Doppler velocity an d temperature 


calcul ation with CRISP were performed following the same technique as in lAntolin et ah 


(I 2 OI 2 I) and we refer the reader to that paper for details. It is important to note that with 
this technique upper limits to the temperatures are provided, since the micro velocities from 
turbulence are ignored. 


For th e follow-up of individual rain features we used the CRisp SPectral EXplorer 
(CRISPEX, Vissers fc Rouppe van der Voortll2012 ') and its auxiliary program TANAT (Times- 
lice ANAlysis Tool). Both are widget-based tools, which enable the easy browsing and anal¬ 
ysis of the image and spectral data, the determination of clump paths, extraction and further 
analysis of space-time diagrams. Fitting piece-wise segments along the clumps trajectories 
in the space-time diagrams we can estimate the projected velocities (in the plane of the sky). 
With this method, the error (standard deviation) for projected velocities is estimated to be 
roughly 5 km s“^. 


3. Temperature evolution 
3.1. SDO/AIA & 55T/CRISP observations 

In Fig. [1] and Movie 1 active region AR11903 is viewed at the East limb. The hgure 
is a composite of three AIA hlter images, combining the 193 hlter (in red), dominated by 
the coronal Fe XII 193.509 A emission at logT = 6.2, the 171 hlter (in green), dominated 
by the Fe IX 171.073 A emission at logT = 5.9, and the 304 hlter (in blue), dominated by 
the transition region He II 304 A emission at logT = 4.9. The movie clearly shows that 
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a significant quantity of the coronal volume in the active region, hlled by loops, exhibits 
continuous intensity var iation. Such dynamical state is strongly suggestive of loops out of 
hydrostatic equilibrium flAschwanden et ahll200ll) . and possibly in a thermal non-equilibrium 
state. This is further supported by the fact that large parts of the coronal volume change 
from red to green and to blue, indicative of cooling. 

During this observation the SST was pointing at the footpoint region of these loops, as 
indicated in Fig. [H where a sunspot is located. The sunspot is roughly located in the middle 
of the FOV. 


3.2. Cooling from EUV to Ho temperatures 

To investigate the thermal non-equilibrium scenario we have chosen a loop (hereafter 
loop 2, which occurs after loop 1, dehned later on) in the large SDO FOV, which exhibits this 
gradual change of intensity from the hot to the cool AIA channels. This loop is highlighted 
in Fig. [T]in dotted lines, it extends ~ 42 Mm above the surface leading to an approximate 
length of 132 Mm for a circular geometry. Its left footpoint is also covered by the FOV of 
the SST. A zoom-in hgure focusing on the loop is shown in Fig. [3l where two snapshots in 
the loop evolution can be seen in the upper two panels, separated by ~ 24 min, while the 
lower panel shows the integrated intensity in different AIA channels along the top part of 
the loop, as well as the size-weighted coronal rain intensity in Ha (dehned later on) close 
to the loop footpoint. These panels show clearly a progression in the AIA channels from 
193 dominated emission to 304 dominated emission. Prior to t ~15min the loop top is only 
visible in 193 but decreasingly so between t = 10 — 17 min. From that point on 171 starts 
dominating clearly until peaking at t = 25.8 min. Only when 171 has already decreased to 
a local minimum plateau do we observe an increase in 304 intensity at the top, which takes 
over from t = 41.0 min and peaks at t=50.6min, close to 25 min after the peak in 171 signal. 
The hrst sign of coronal rain in Ha in this loop is already detected prior to the 304 peak. 

While AIA 304 has a single-peaked light curve, AIA 193 and 171 show much more 
variability on the scale of a few to 10-15 minutes. The variations in Ha are much more 
extreme and on a shorter timescale than for the AIA channels, but this is to be expected 
considering the difference in cadence and the fact that the AIA intensities are integrated over 
a large area, while the Ha intensities stem from single rain clump measurements. A large 
quantity of Ha clumps is observed falling across the SST FOV with basically two slightly 
different trajectories but converging to similar locations, apparently within the umbra of the 
sunspot. In Fig. Ulan example of a large clump is shown in an Ha Doppler image (±0.76 A, 
left panel) and in AIA hlters. The track of the clump, shown in dotted lines coincides with 
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a dark loop structure of larger width in the EUV passbands. This dark structure matches 
with a bright loop structure above the limb (hereafter loop 1), adjacent (to the right) to 
loop 2. A composite image of the 5'5'r/CRISP FOV (Doppler image Ha ± 0.76 A in red) 
together with AIA 171 (in green) and AIA 304 (in blue) is shown in Fig. El where the same 
rain event as in Fig. |4]is shown. 


In order to identify all the clumps, for instance those belonging to loops 1 and 2, a 
semi-automatic routine was written detecting all rain-like structures crossing a specihc path 
in the FOV (see Section [23]) • This routine calculates the intensity, Doppler velocity, the 
width and the length for each clump. We select two cuts. Cl and C2, roughly transverse 
to loops 1 and 2, respectively, plotted in Fig. El crossed by most of the clumps. Movies 3 
and 4 show the detected clumps as they cross Cl and C2, respectively (in red the FWHM 
of the £t to the clumps’ widths). These cuts can be seen to be roughly perpendicular to the 
trajectory of the clumps. The algorithm conhdently detects 46 and 95 clumps for Cl and 
C2, respectively (as shown in the movie, a myriad of smaller clumps exist, but these are too 
faint and cannot be conhdently detected by the algorithm). Despite the apparent crossing of 
these loops (evidenced by the intersection of Cl and C2), these are actually well separated 
along the LOS (the algorithm further ensures that no clumps are common to both Cl and 
C2). This is evidenced by the different Doppler velocities between clumps belonging to 
separate loops. Indeed, clumps belonging to a specihc loop will share a common trajectory 
and therefore a similar Doppler velocity in general. This correlation is also obtained in 
3D num erical simulations of thermal instability reproducing coronal rain and prominence 
threads fiLnna et ah! l2012l) . Figure El shows a histogram of Doppler velocities and derived 
temperatures for all detected clumps. Positive and negative values correspond, respectively, 
to redshifts and blueshifts. We can see that the Doppler velocities concentrate in 3 diherent 
groups, one peaking at ~ 30 km s“^, another at ~ —10 km s“^ and the last one at ~ 
—40 km s“^. The hrst and second groups correspond to the clumps crossing Cl and C2 
respectively, while the third corresponds to the Doppler velocities of the clumps in loop 2 
at chromospheric heights (cut F2, defined in Section 14.41) . The projected velocities of the 
clumps crossing Cl and C2 are found to be similar, with an average of 110 km s“^. This 
ensures that the angles of fall for each family of clumps is diherent, and similar within each 
family (as indicated by the formulae for the falling angle calculated in Paper 1). We conclude 
therefore that the Cl and C2 clumps correspond to 2 diherent loops. 


The clumps crossing C2 are found to follow a dark structure in the EUV passbands 
similar to that plotted in Fig. [3 This dark EUV structure appears bright above the limb 
and matches with loop 2. The geometry of the loop seen in Fig. Eland the direction of the 
how suggest that the material is coming towards the observer, matching with the negative 
Doppler velocities found in Fig. |6]for the C2 clumps. We can therefore conhrm that the C2 
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clumps belong to loop 2. 


As shown by lAnzer fc Heinzell (120051) the Ha intensity from hlaments is strongly cor¬ 
related to darkening in EUV pass bands, and is due mainly to continuum absorption from 
neutral hydrogen, neutral helium and singly ionised helium (and to a lesser extent to volume 
blocking, especially in the case of small size structures as coronal rain clumps). Therefore, 
both the volume and the composition of the emitting structure are important factors that 
can be correlated to the detected intensity variations in the EUV passbands. In the lower 
panel of Fig. Elwe thus plot as a black histogram the Ha intensity multiplied by the width 
squared (assuming axi-symmetry for the clump) for the clumps crossing C2. The bulk of the 
distribution spans over 25 minutes and the peak is shifted by ~ 20 min from the 304 peak. 
However, the Ha detection is made at the footpoints, while the AIA peaks are hrst observed 
at the apex. The clumps are observed to fall with average total speeds of 110 km s“^ and 
accelerations around 0.1 km s“^. Assuming no initial velocity, constant acceleration and a 
circular loop, the travel time of the cooling plasma is 20 ± 10 min, which locates in time the 
origin of the clumps close to the 304 peak. This result supports the common origin for the 
Ha clumps with the rain observed in 304 in loop 2. 


The temperature histogram in Fig. |6] shows that all the detected clumps have very 
cool temperatures, mostly below 10000 K, with a peak at 2000 — 5000 K. As explained in 
Section 12.31 these temperatures correspond to upper limits. While we postpone a discussion 
on these low temperature values to the pertinent section, these results clearly indicate that 
loops 1 and 2 undergo full catastrophic cooling to chromospheric temperatures. 


3.3. EUV variation associated with Ha rain emission 

It is interesting to compare the EUV and Ha intensity variations across the Cl and 
C2 segments. This is shown in Fig. [71 For each time step we average the intensity in each 
channel over a small region centred around each cut. For Ha, as in Fig. [HI we plot the 
average of the intensity times the square of the width for each clump (see hgure caption for 
more details). Also, for the times at which a clump is detected the intensity in each EUV 
channel is calculated over the location of the clump only, in order to discern better a possible 
effect on the EUV passband for the passing clump. For the clumps crossing Cl we can see 
that most of the time groups of clumps close in time produce a decrease in all of the EUV 
intensity channels (at times t ~ 0,10, 27,40 min), except 304, for which an intensity increase 
is generally detected. For C2 a clear correlation is also detected. Quasi-periodic fluctuations 
in the EUV passbands observed in the time range t=60 — 90 min seem to match the Ha 
intensity fluctuation in the same time range. Contrary to the previous case here a maximum 





12 


in the Ha intensity seems to be positively correlated with all EUV passbands. Indeed, each 
peak of the former seems to correspond to a local peak in the latter. 


3.4. Hinode/SOT, IRIS/SiY Sz SDO /observations 

We now tnrn to the second set of observations involving SOT, SJI and AIA. A look at 
Fig. [2] and Movie 2 immediately indicates the large presence of cool material above the active 
region. In red, green and bine we show Ca II H emission from SOT, Si IV 1400 emission from 
/i?/S'/SJI and Fe IX 171 emission from AIA, respectively. We can see the presence of two 
kinds of cool strnctnres above the spicular layer. Prominence material is observed flowing 
relatively slowly ronghly horizontally, while coronal rain is observed falling relatively faster 
along loops bright in AIA 171 (other cool static prominence-like strnctnres are also observed, 
e.g. at {x,y) ~ (—160,970)). Here we will concentrate on a loop exhibiting coronal rain, 
with one footpoint at {x,y) ~ (—120,930) and the other apparently at {x,y) ~ (—150,955), 
expanding in the corona np to a height of ronghly 40" (leading to a length of ronghly 180 
Mm assnming circular geometry). The rain appears from 23:12 UT to 23:30 UT (the end of 
the co-observation sequence) and is observed to initiate at the apex of the loop and to fall 
towards the farther footpoint. 

A zoom-in onto the loop of interest is shown in various wavelengths in Fig.[H]and Movie 7. 
Contrary to the hrst dataset, here we are able to observe the chromospheric emission from 
the rain well above the spicular layer, but we are not able to follow it once it goes below this 
layer. 


3.5. Cooling through transition region temperatures 

In order to follow closely in time the evolution of temperature in each wavelength we 
plot in the upper panel of Fig. [9] the integrated emission for each hlter for coronal rain 
pixels belonging to the previously indicated loop. As seen in Movie 2, half-way through the 
coronal rain event a prominence that appears to be in the background crosses the upper 
portion of the loop. In order to minimally reduce such LOS projection effects the selection 
of clump paths and the times and locations of rain along such paths is done manually and 
is done ensuring that the prominence is not along the LOS. In the hgure we can see that 
the hrst emission of rain appears in the /i?/S'/SJI images, before showing up in the Ca II H 
hlter. The SJI 1400 hlter comes hrst (at to ~ 4:4.5 min), followed by a peak in SJI 1330 
(f = to + 2.5 min), then SJI 2796 (f = to + 5.5 min) and a small peak in Ca II H roughly at 
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the same time. The light curves generally show short-timescale variability. While all the SJl 
channels have a hrst large peak followed by gradually smaller peaks, the Ca II H emission 
starts with a small peak and increases over the next 5 minutes, with the highest peak at 
t = to+ 8 min. Interestingly, the time lag between the hrst or highest peak of each diagnostic 
seems to increase as the rain becomes visible in each: the C II peak follows the Si IV peak 
after ~2.5 min. Mg II k follows C II after ~3 min, Ca II H follows Mg II k after ~ 3 — 3.5 min. 
Although the dynamic range differs, part of the increases and decreases are visible accross 
several diagnostics. For instance, the initial C II peak around 46.5 min seems to coincide 
roughly with a secondary Si IV peak, the initial Mg II k increase and peak around 49 min 
is mirrored in minor by Ca II H, while AlA 304 and I7I seem to follow each other closely 
throughout. 

The behaviour of the AlA hlters, and especially 171, is strikingly different. Movie 2 
shows that most of the loop is bright and stays bright in I7I, from a hrst peak about 10 min 
before hrst appearance of coronal rain in SJl 1400. No EUV darkening co-spatial to the 
rain is observed in this case, as was the case for dataset 1. Both AlA 171 and 304 increase 
slightly in intensity during the rain event, with the 304 emission of the rain starting at 
t = to + 3.5 min, roughly at the same time as the Mg II k emission. This is more clearly 
seen in Movie 7. This behaviour seems contrary to that observed in the loops in the hrst 
dataset with CRISP and AlA, in which 171 emission at the loop apex peaks either prior or 
following a rain event and is at a minimum during the rain event, and in which 304 peaks 
after 171 and roughly at the same time or prior to Ha (cf. Fig. [3]). Multiple small peaks 
are also observed in the AlA 171 emission, in a quasi periodic fashion separated by a few 
minutes, similar to the behaviour observed in dataset 1, shown in Fig. [3 

Another useful way to visualise the evolution of emission is through the cumulative light 
curves, plotted in the lower panel of Fig. O In this hgure we can see the same progression 
throughout the hlters, with the variability rehected now in the slope of each light curve. 
Accordingly, the behaviour between the SJl hlters is very similar, characterised by two 
diherent slopes: an initial fast rise with a steep slope (until t 50 min) followed by a more 
gradual rise. The fast rise can also be seen in movies 2 and 7: the rain brightens as soon 
as it appears close to the apex of the loop, especially in SJl 1400 and 1330. This two step 
behaviour suggests two diherent cooling phases. The Ca II H emission appears mostly during 
the second cooling phase, it has roughly the same slope throughout the evolution, rehecting 
a similar rise and decay of the intensity in the upper panel. Accordingly, movies 2 and 7 
show a more gradual emission in Ca II H. The behaviour between the AIA 304 and 171 
hlters is very similar and constant on average, leading to a roughly constant slope in the 
cumulative plot of Fig. [9l 
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Figs. [HI [9] and Movie 7 further show that the emission in all hlters co-exists in the loop 
during most of the catastrophic cooling event, but especially during the second phase of slow 
cooling. This strongly suggests a multi-temperature structure for the thermally unstable 
loop. This is further supported by the increase of the AIA 171 intensity along the loop 
during the rain event. This will be further discussed in Section [521 

The apparent mismatch with the AIA 304 channel, this one occurring 3.5 min after hrst 
occurrence in SJI 1400, could be due to the low sensitivity and lack of resolution of AIA 
with respect to the other instruments. A small clump will not become bright enough above 
the noise level unless it either increases in intensity or becomes larger. In Movie 7 we can see 
that AIA 304 emission is barely visible above the noise level, supporting the interpretation. 


4. Morphology 

4.1. Rain showers with cool cores and associated continuous EUV darkenings 

Figure 12] and Movie 1 show that the downward flow along many of the loops in the active 
region seen in AIA 304 has a continuous character with little clear substructure, especially 
in the direction of motion. This is also observed in the lower part of the loops shown in 
Fig. m where the downflow appears dark in absorption. However, as revealed by movies 3 
and 4, the falling Ha clumps observed with CRISP come in various sizes and intensities along 
the same paths. Fig. [2] shows an example of a clump in Ha in emission above the spicular 
layer, and extending onto the disc towards the umbral region of the sunspot (marked by a 
dotted line in the hgure). The on-disc part of the elongated rain clump is observed as a 
thin absorption feature and is adjacent to other elongated rain clumps (also seen in Fig. |4|). 
Surrounding the Ha emission, EUV darkening can be observed all along the clump trajectory 
in the lower part of the loop. This low EUV emission is more persistent, as shown in Fig. [T] 
thus matching the continuous (as opposed to clumpy) character observed in AIA 304 in 
the larger FOV of Fig. [1] However, as shown in Fig. [7| and described in Section 13.31 a 
close look at this low emission reveals small EUV variations, which appear correlated to 
the Ha intensities. This occurs especially when several Ha clumps appear close in time. 
Such groups composed of a large number of small Ha clumps were denoted in Paper 1 as 
‘showers’. As predicted in that work here we show that such showers, and also a small 
minority of large and dense enough Ha clumps, are associated with absorption features in 
EUV (and corresponding intensity variations within the dark structure). Figure |4| shows 
that these absorption features are signihcantly wider than the Ha emission (on the order of 
a few times the Ha width), suggesting a wide transition from cooler to hotter temperatures. 
Also, a comparison between movies 1 and 3 (or 4) indicate that the EUV absorption features 
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are significantly longer than the Ha clnmps (except for the few very long clnmps that extend 
beyond the FOV, for which we cannot conhrm). These results may suggest at hrst a scenario 
in which cool chromospheric cores are surrounded by warmer diffuse shells. However, the 
large difference in spatial resolution between CRISP and AIA may be the main factor behind 
this picture. 


4.2. Co-spatiality of emission and substructure 

The previous scenario is further supported by dataset 2. As shown by Fig. [S] and 
Movie 7, the downflow in AIA 304 appears mostly structure-less, although with a slightly 
brighter head (which also appears bright in the other chromospheric and transition region 
hlters). On the other hand, substructure can clearly be seen not only in the chromospheric 
hlters (Ca II H and SJI 2796) but also in the transition region hlters (SJI 1330 and SJI 1400). 
Close inspection further indicates a very similar rain structure suggesting a high degree of 
co-spatial emission in both chromospheric and transition region hlters. This can more clearly 
be seen in Fig. [TOl which shows a composite image of Ca II H (in red), SJI 1330 (in blue) 
and SJI 1400 (in green). This result suggests that the lack of structure in AIA 304 may not 
be due mainly to a diherence in temperature but rather to a lack of spatial resolution. 

To analyse more closely the degree of co-spatiality we zoom-in into the region and 
follow a small shower of clumps (the head of the rain, which appears bright in Fig. |8]). 
This is shown in Fig. [TTl where the same shower is shown at three diherent times along 
its fall, separated by ~ 40 s from each other. The shower’s length becomes shorter as 
it falls (passing from 8" to 4"). Its width on the other hand remains mostly unchanged. 
However, the width of the shower (and the amount of substructure) depends strongly on 
the wavelength. While in AIA 304 the shower appears as a single entity, the presence of 
substructure is obvious in Ca II H. The shower can be seen to be composed of two strands, 
each one presenting inhomogeneities along their lengths. The substructure is far less evident 
in SJI 2796, probably due to the signihcant increase of opacity in this line. However, it 
resurfaces in SJI 1330 and SJI 1400 with similar shapes, although less clear due to the lower 
resolution. This picture not only proves a multi-temperature scenario for coronal rain but 
also proves the existence of substructure in thermally unstable loops in chromospheric to at 
least transition region temperatures. A look at AIA 171 (Movie 2) further shows that coronal 
temperatures are present throughout the loop, especially in the wake of the rain head. 
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4.3. Density for a thick clump 


Through the change in AIA emission produced by EUV absorption from hydrogen and 
helium in the c l umps in dataset 1 we can estimate their densities following the technique by 
Landi fc Realel (120131) . The application of this technique is based on certain assumptions. 
Mainly, the plasma needs to be in ionisation equilibrium, the produced absorption needs to 
be homogeneous along the LOS, and the emission in the corresponding wavelengths of the 
absorption needs to be minimal. While these may be questionable in the present case (a 
catastrophically cooling plasma could be out of ionisation equilibrium if the cooling timescale 
is shorter than the ionisation and recombination timescales; also, as shown previously, we 
have a strong degree of density inhomogeneity within the rain) the results of the applied 
technique can nonetheless provide important density constraints and clues on the state of 
the plasma. Here we do this for a large clump producing enhanced EUV darkening, similar 
to that shown in Figs. |l]and[5l but we postpone a full investigat i on on densities to future 
work. In Fig. [T2]we plot the L function, defined in lLandi fc Realel (120131. formula 9), which, 
at temperature of absorption is equal to the hydrogen column density. Our case is similar 
to their special case 2 (Section 3.2 in their paper). The background emission is measured 
at the same locations where the EUV darkenings are observed but at times where no Ha 
clumps are observed (we take a few time steps prior and after the passage of the clump, 
and then take the average of the two). The EUV intensities are taken over the trajectory 
of the clump, seen as darkened streaks in an x — t plot. In the calculation of L we further 
assume a helium abundance of 5%. In Fig. [12] we can see that the region of intersection of 
all curves is concentrated in logTabs = 4.4 — 4.6 and logL = 18.1 — 18.7 cm“^ (where Tabs 
is the temperature of maximum absorption). Assuming axi-symmetry along the direction of 
propagation the LOS depth across the clump is then roughly equal to the clump’s width. 
The clump’s width in the EUV filters is 2.8", although it cannot be determined with precision 
due to the poor spatial resolution of AIA. The clump in Ha has a maximum width of 0.6" 
(with an average of 0.4"). From the results of dataset 2, it is likely that the larger widths 
in EUV are mostly an effect of the lack of resolution (which would smear out the width a 
factor proportional to the AIA’s PSF). We therefore assume an EUV width of 1.5 — 2 times 
larger than the Ha width, leading to a width of roughly 700 km. This leads to an electron 
density of 1.8 — 7.1 x 10^° cm“^, where the scatter is most likely due to the presence of 
inhomogeneities within the clump. This is supported by the fact that the spectral width of 
the clump in Ha indicates a temperature of 5500 ± 500 K, suggesting a multi-temperature 
structure, and possibly a higher density. Assuming a constant pressure within the clump the 
core density can then be up to 2.5 x 10^^ cm“^. 
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4.4. Clumpy vs. continuous 

As the clumps fall further down towards the sunspot in dataset 1 they become absorption 
features in the wing of Ha while they gradually become bright towards line centre (see 
Figs. |l]and[5]). These distinctive features allows us to trace them down to chromospheric 
levels, just before impact. Similarly as previously done for coronal heights we dehne two 
cuts at chromospheric heights roughly perpendicular to the clumps trajectories, FI and F2, 
corresponding to loops 1 and 2, respectively. Both C & F pairs are plotted in dashed curves 
in Fig. |5l As movies 3 and 4 did for Cl and C2, movies 5 and 6 show the detection of rain 
clumps across FI and F2, respectively (in red the FWHM of the ht to the clump’s width). 
These movies indicate that the Cl clumps fall either into the umbra or towards the closer 
edge of the sunspot (projection effects do not allow us to hrmly pinpoint the location of the 
fall) and the C2 clumps seem to fall mainly towards the far edge of the sunspot. Tracing of 
coronal rain clumps to such low heights has never been achieved before and is possible here 
thanks to the high falling speeds of the clumps (which produces the dark absorption features 
in the wing of Ha) and to the high spatial and temporal resolution of CRISP. Indeed, the 
total velocities (taking into account projected and Doppler velocities) vary between 80 km s“^ 
and 120 km s“^ for both loops. 

Due to the increase of projection effects at low heights, the high speed of the clumps 
and the seeing effects (which not always provides clear subsequent images), it is not possible 
in most cases to specify the start and end for a clump at low heights, and therefore provide a 
robust determination of clump lengths. It is therefore possible that two subsequent detections 
at the same location correspond to the same clump. Also, it is highly likely that many of the 
clumps detected at heights C are detected again at heights F. However, since we are more 
interested in regarding the rain as a flow rather than a discrete set of condensations (especially 
because their morphology can change dramatically as they fall), for our statistical analysis 
we consider that all clumps are different (and for the lengths we give an approximate value 
determined visually). Since there are many clumps that fail to be detected by the algorithm 
(due to the previously stated reasons), we consider that this procedure leads to a better 
approximation to the real rain population. 

A striking feature shown by movies 5 and 6 is that at such low chromospheric heights 
the clumps appear signihcantly less clumpy (their lengths increase) and the downward flow 
becomes continuous and rather persistent. In order to quantify this effect better we plot 
in Fig. [13] a histogram of clump detection weighted with their respective widths squared at 
the two previously dehned heights for loops 1 and 2. Not only the clumps in the corona are 
detected at later times in the chromosphere, but we also detect many more rain crossings 
at those lower heights, indicating an increase in the rain flow with decreasing height (46 
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and 145 clumps for Cl and FI, respectively; 95 and 193 clumps for C2 and F2, respec¬ 
tively). In Section [53] we discuss possible reasons for this increase in clump occurrence at 
low chromospheric heights. 


4.5. Strand-like structure 


As shown in Section 14.21 the presence of substructure in coronal rain at high resolution 
is not only clear at chromospheric temperatures such as those of Ca II H and SJI 2796 but 
also at higher temperatures, such as those of SJI 1330 and SJI 1400. As shown in Fig. [TTl 
as the rain falls substructure is observed to be organised along strands of various lengths 
but rather similar widths. This is also observed in Movie 7 for the rest of the rain in this 
event and in dataset 1 at even higher resolution. Indeed, at the highest resolution achieved 
in this work (in Ha with CRISP) similar strand-like structure appears, as shown in Fig. |4| 
and movies 3 — 6. Due to the very dynamical nature of coronal rain it is not easy to discern 
the strand-like structure traced in time by the rain. For this, we plot in Fig. (TT] the time 
variance in a Doppler image in Ha ± 0.6 A over a time period of 22 min for loop 1 and 
its surroundings (where the time variance is dehned as the sum of the squared average- 
subtracted intensity for each pixel). Thanks to the emission features of the rain in the wing 
of Ha both above the spicular layer and at chromospheric levels, the clumps produce bright 
traces in the variance image. These traces appear as various strands at both of these heights 
(in the hgure, these are visible from 50"to 55" and 28" to 35" in the vertical axis). In the 
hgure, 6 to 8 strands with widths around 0.2" — 0.4" can be discerned at both heights, 
although a larger number may exist due to projection effects. It is worth noting that the 
22 min period over which the variance is taken is a signihcantly long time for coronal rain. 
Indeed, the number of clumps detected over this time period is on the order of 60 — 70. Also, 
as calculated previously in Section 1331 the travel time for rain in these loops is in the order of 
20 ±10 min. Lastly, numerical simulations indic ate that the catastrophic cooling part of limit 


cycle s is generally short, under an hour or so flAntolin et al.ll20I0l: iMendoza-Briceho et ah 


20051: ISusino et ah! l20I0h (although this depends on various parameters of the heating). 


Therefore a signihcant part of the catastrophic cooling for this loop may be captured in the 
time variance, strongly suggesting that the clumps follow well-dehned fundamental magnetic 
structures in the corona. 


Interestingly, apart from the prominent clumps as those in Fig. [TTl fhe surrounding 
diffuse atmosphere in dataset I, dimly bright in Ha also appears structured in the same way. 
This is especially clear at times of large clumps. An example of this is shown in Fig. [T5l 
where a snapshot of a falling clump can be seen in Ha ± 1.1 A. A bright clump with a 
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long tail can be seen and, particularly, multiple strand-like structures next to it, as ripples 
with similar widths down to 0.15" or so, uniformly distributed within 4" of the clump. This 
structure can still be discerned (although with decreasing intensity) when shifting the wave¬ 
length position over a range of ±0.5 A around 1.1 A, ruling out possible instrumental effects 
from the narrowband CRISP wavelength channels. In order to show the strand structure 
better we make transverse cuts across the clump at two different locations, one at the head 
and one at the tail of the clump. The H« intensity along these cuts is shown in the right 
panels of the hgure. Next to the large peak corresponding to the bright clump, around 8 
ripples can be distinguished extending significantly along the direction of propagation. This 
structure is highly reminiscent of the MHD thermal mo de (also know n as the entropy mode 
in the absence of thermal conductio n), first predicted bylFieldl (1965) and later developed b y 
van der Linden fc Goossens (1991b a); Goedbloed fc PoedtsI (12004 '): Murawski et ah (12011 1. 


Indeed, the spatial distribution of such mode is that of a main elongated dense clump with 
multiple smaller clumps side by side of similar widths. The generation of this wave is guar¬ 
anteed by the small but non-zero perpendicular thermal conduction in the corona. Although 
static in an ideal scenario, this wave is expected to move together with the flow. In this sce¬ 
nario, the observed strand-like structure at the smallest detected scales could be the result 
of MHD thermal modes produced from thermal instability. 


4.6. Sizes 

We now turn to the statistical determination of sizes. The methods used in the widths 
and lengths calculations are explained in Section [2731 The routines used for this detection are 
semi-automatic. The detection of a clump is rather strict, and is based on several conditions 
which rule out most of the fainter and smaller clumps at diffraction limit resolution. For 
instance, the smaller widths that can be discerned in the cross prohles in Fig. [15] are mostly 
left out by the algorithm, in favour of precision and reduction of errors. 

In Fig. [in] histograms displaying the relative frequency (for each wavelength) for the 
obtained widths and lengths are shown, combining datasets 1 and 2, and therefore all the 
instruments used in this study (note that the bin widths of the distributions have been 
set equal to the resolution of each instrument). Due to the large number we restrict the 
histograms to the more interesting small scales part {upper panels) and separate between the 
high (CRISP, SOT, SJI) and low resolution instruments (AIA in the lower panels). The hgure 
clearly shows a signihcant difference between widths and lengths at all wavelengths. The 
widths appear nicely distributed, with increasingly higher occurrence frequency at smaller 
scales for each wavelength and a tail at longer scales. Such distribution shape seems therefore 
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independent of temperatnre, with an apparent peak for all clnmps between 0.2" and 1.1". 
Ha and Ca II H clump distributions have, respectively, narrow and broad peaks at 0.25" and 
0.6". Blobs observed with IRIS in SJI 2796, SJI 1330 and SJI 1400 have similar distributions 
peaking at 0.6" — 0.8". We note however that SJI 2796 presents a broader distribution for 
the peak, between 0.6" and 1", probably due to the higher opacity compared to the optically 
thinner C II and Si IV lines. For AIA we note that while widths in 171 and 193 peak around 
1.8" — 3", 304 peaks around 0.8" — 1.7". 

On the other hand, the lengths show a more random distribution, although with high 
clustering at small numbers, between 0.5" and 1.5". It is worth noting that the lengths 
distribution spans up to 30" or so (for instance, the very long clump of Fig. 0]). For lengths 
longer than 5" (not shown here) the distribution becomes sparse. 

In Table 0] average values for each wavelength and dataset are shown, together with the 
total number of clumps detected in the FOV per minute as an indication for the number of 
clumps detected in each dataset. We choose to calculate this number since a clump is a not 
well dehned plasma structure over time. As shown in the present movies (and in Fig. [TT]) 
and in Paper 1, clumps can dramatically change shape, appear and disappear along their 
fall. This number needs to be read with caution: the area over which rain is detected in 
dataset 1 is roughly half that in dataset 2 (therefore, dividing the number over area would 
lead to a similar number between SST and Hinode). The continuous flow in 304 (as opposed 
to clumpy) also makes it difficult to identify a clump. For AIA only clear cut examples of 
clumps were selected, which explains the low number per minute. 


The widths presented here in Ha are similar to previously measured widths in the same 


wave length with CRISP in large statistical datasets (lAntolin et al.ll2012l:lAntohn fc Rouppe van der Voort 


20121 ). but present a relatively higher number of clumps at lower spatial scales. This is mainly 


due to the fact that the CRISP instrument gained from improved spatial sampling after an 
upgrade perfo r med in 2009. The widths of Ca II H clumps are similar to those reported in 
Antolin et al.l (120101) . In the top left panel of Fig. [16] the tendency of number increase for 


clumps with smaller widths suggests a tip-of-the-iceberg scenario. 


As expressed previously, the lengths appear more randomly distributed, especially in 
Ha, and therefore may be more dataset dependent. For instance, extremely long clumps are 
only observed in dataset 1, and have not been previously obser ved in other datasets in Ha 
with CRISP, except in the recent work by IScullion et al.l (120141) in the context of flares, i.e. 
stronger footpoint heating sources. 
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Table 1. Average values for both datasets 


Instrument 

& hlter 

Dataset 

Width 

[arcsec] 

Length 

[arcsec] 

Blobs 
/ min 

CRISP - Ha 

1 

0.30 ± 0.09 

3.85 ± 5.77 

5.0 

SOT - Ca II H 

2 

0.61 ± 0.17 

1.05 ± 0.51 

9.7 

SJI - 2796 

2 

0.80 ± 0.20 

2.49 ± 1.97 

2.6 

SJI - 1330 

2 

0.70 ±0.15 

1.75 ± 1.35 

3.7 

SJI - 1400 

2 

0.81 ±0.17 

1.76 ± 1.26 

4.6 

AIA - 304 

2 

1.56 ±0.15 

10.8 ± 6.5 

2.2 

AIA - 304 

1 

1.40 ±0.70 

12.7 ± 9.2 

1.2 

AIA - 171 

1 

2.45 ±0.81 

5.5 ± 5.2 

1.2 

AIA- 193 

1 

2.69 ±1.07 

7.4 ± 4.6 

1.2 


Note. — Averages are calculated over all detected blobs in 
both datasets. Row order is set according to increasing tempera¬ 
ture. For dataset 1 the measurements are limited only to the rain 
detected in loops 1 and 2. For dataset 2 the measurements are 
limited to the loop shown in Fig. [HI While each measurement is 
itself an average value over a large quantity (see section 12.3p and 
therefore has its own standard deviation, the standard deviation 
shown here for widths (lengths) is simply taken over the set of 
blob widths (lengths). The last column provides the number of 
blobs detected per minute for each dataset, which implies depen¬ 
dence over the FOV over which this detection is carried. For the 
SST, the detection is done over a small part of its FOV (roughly 
5"xl7"), which is much smaller than the FOV for the SDO data 
of the same dataset (see Fig. |3]). For Hinode the projected area 
over which the the rain is detected is roughly twice that of SST. 
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5. Discussion 

5.1. EUV intensity variations as a signature of catastrophic cooling 

Coronal loops in active regions often present intensity variations, which are regularly 
associated to states out of hydrostatic equilibrium such as states of thermal non-equilibrium 
flAschwanden et ahll200ll: iRealel 120101) . A state of thermal non-equilibrium implies the pres¬ 
ence of thermal instability in coronal loops. This instability can be thought as the most 
general state of the plasma (meaning that thermal equilibrium is extremely unlikely). The 
relevant question here is what are the timescales of this instability in coronal loops? Accord¬ 
ingly, this instability can have a complete character, in the sense that runaway cooling due 
to the shape of the optically thin loss function may occur locally in coronal loops and con¬ 
tinue uninterrupted down to chromospheric temperatures, until heating from local sources 
overcome the losses due to its own radiation. The completeness of thermal instability is, 
however, not guaranteed an d depends st r ongly on the efficiency of the heating mechanisms 
and on thermal conduction. iMikic et ahl (120131) have shown that even changes in geometry 
(non-uniformity in area cross-section and asymmetry in the heating) play an important role 
in determining the aspect of this instability. Furthermore, the spatial scales involved in the 
cooling processes (for instance, the dumpiness and strand-like structure of the rain) are di¬ 
rectly related to the properties of the plasma and to the counter-acting heating mechanisms, 
and may therefore reveal essential characteristics of the heating processes. 

Despite the close link between thermal instability and coronal heating and therefore the 
importance of thermal instability, we do not know how ubiquitous this instability is in the 
corona and the details of its character (for instance, whether it is complete or not). This is 
largely due to the high complexity of the thermal instability in strongly anisotropic plasmas 
such as the solar coro na. Recent observationa l studies over active regions indicate a general 
tendency for cooling fjViall fc Klimchukll2012l) . Although it is clear that not all the loops in 
an active region are in a thermal non-equilibrium state and undergo catastrophic cooling, it 
is important to estimate what fraction of loops are. 

The traditional picture of thermal instability implies that coronal loops in a thermal 
non-equilibrium state become progressively cooler, suggesting a time delay between hlters 
detecting emission at different temperatures. Such cooling progression throughout transition 
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since all of these multi-wavelength observations do not include chromospheric ranges. In this 
































































23 


work we make the first steps in this direction by linking commonly found intensity variations 
in EUV hlters of SDO/AlA to the observational imprint of catastrophic cooling, i.e. coronal 
rain observed in chromospheric diagnostics. 


We observe complete thermal instability with temperatures down to the chromospheric 
range. In Section[3]we show that EUV darkening is found to be associated to flows in AIA 304 
and presents a continuous and rather persistent character. Such AIA 304 flows are observed 
along many of the loops within the active region (dataset 1), part of which are captured 
towards the footpoints in Ha by CRISP. The flows in AIA 304 appear strongly correlated in 
space to the appearance of cool chromospheric material in the loops observed in Ha, which 
occurs in a more intermittent way. This suggests that while the entire loop may be thermally 
unstable and cooling, plasmas at transition region temperatures are more widespread and 
common within the loop than plasmas at chromospheric temperatures. This is expected to 
some extent due to the tendency of the plasma to keep a uniform pressure distribution along 
the held (and the fact that thermal instability entails a local loss of pressure leading to small 
condensations). 


Small EUV intensity perturbations are further strongly correlated in time with the 
presence of Ha rain clumps. Although intermittent, these clumps can come in series of 
events and generate quasi-periodic EUV intensity variations on the order of a few minutes. 
In most cases, these variations are found to be associated to groups of chromospheric clumps 
which we have previously termed ‘showers’ (Paper 1), and span over a few arcseconds in 
the transverse direction. Correlation between EUV darkeni n g and Ha has been e xtensively 
studied in the case of prominences flHeinzel fc Anzerll2006l: iLabrosse et al.N20101) . For the 
EUV wavelengths of interest here, this correlation is mainly due to continuum absorption 
from neutral hydrogen, neutral helium and singly ionised helium. 


5.2. Multi-thermal character 

The progressive cooling picture from thermal instability also suggests a difference in 
emission with height. In this picture, emission in Ha and at cooler wavelengths is preferen¬ 
tially detected at low coronal heights. In the middle panel of Fig. [T7]we plot the averages 
of height versus temperature for all detected clumps in datasets 1 and 2. The temperatures 
taken for the SJI and SOT hlters are those of maximum formation for the dominant ions of 
each hlter. Those for CRISP correspond to the measured temperatures, shown in Fig. |6l Al¬ 
though not strongly correlated, the plot suggests a slight temperature decrease with height. 
Particularly, the very low temperatures below 5000 K, even down to 2000 K or so, for the Ha 
clumps detected with CRISP at the lower part of the loops provides support to the cooling 
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picture. 

Through co-observation between IRIS and Hinode/SOT we are able to follow for the 
hrst time the rapid progression in the runaway cooling in thermally unstable loops. A two 
step progression is found in which fast sudden appearance at transition region and upper 
chromospheric temperatures is followed by a more gradual appearance at lower chromo¬ 
spheric temperatures. This two step scenario agrees with the shape of the optically thin 
loss function, which predicts runaway cooling down to 10^ K or so, below which material 
becomes progressively optically thick, and therefore radiation has a lower escape probability 
(leading to less efficient radiative cooling). 

Despite the progressive cooling from transition temperatures to chromospheric temper¬ 
atures found in dataset 2 the AIA 171 intensity along the loop is found to increase about 
10 minutes prior to the rain event, and remains bright and even further increases at the 
end, during the rain event. Such behaviour strongly suggests the presence of multi-thermal 
plasma in the loop, from chromospheric to coronal temperatures. Furthermore, it implies 
that the progressive cooling pictnre down from coronal temperatnres to chromospheric tem¬ 
peratures is not always reflected in the evolution of EUV light curves. Hot material can 
co-exist with catastrophically cooling plasma. Therefore a hot to cool systematic sequence 
in EUV channels is not an observable requirement for a loop in a thermal non-equilibrium 
state. 


The obtained cooling down to chromospheric temperatures is very fast, achieved in a 
timescale that varies from tens of minutes (dataset 1) down to minutes (dataset 2). Due 
partly to this reason, in both analysed datasets, and especially in the one analysed with 
SJI and SOT, thanks to the high spatial and temporal resolution we find a high degree 
of co-spatiality in the multi-wavelength emission over the entire span of the event. This 
strongly supports the presence of multi-thermal plasma in these loops. Furthermore, the 
multi-thermal character is accompanied by strong density inhomogeneity, which suggests 
a complex thermodynamic evo lution within the sam e loop. The multi-thermal picture is 
also supported by the results of IScullion et ahl (120141) . which show Ha emission with CRISP 
co-located with EUV AIA 171 in post-flare loops. 


The large range of temperatures we detect imply that coronal rain, as a resnlt of catas¬ 
trophic cooling, is not solely a chromospheric phenomenon bnt a transition region phe¬ 
nomenon as well. Actually, since it is material cooling from coronal temperatures and, as 
discussed previously, it is more pervasive in transition region temperatures, it is mainly a 
transition region phenomenon which spans down to chromospheric temperatures if the ther¬ 
mal instability is complete. This result matches well the picture that is being perceived with 
increasing frequency in SJI and AIA synoptic observations, namely, that basically every ac- 
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tive region shows rain-like phenomena in AIA 304 and SJI 1400. A proper statistical analysis 
of this picture is left as future work. 


For dataset 1, the EUV emission (darkening or brightening) in AIA 304 is found to be 
wider than the chromospheric emission in Ha or Ca II H suggesting a wide transition from 
chromospheric to transition region temperatures. However, the highly co-spatial emission 
obtained from SJI and SOT in dataset 2 suggests that the wide transition in dataset 1 is 
mainly due to the lack of resolution in AIA. The observed thin transition in dataset 2 must 
exist on spatial scales below the resolution of IRIS, i.e. 0.33", which is also supported by the 
similar width distributions obtained with SOT and the SJI hiters (cf. Fig. [16]). Indeed, when 
closely comparing the substructure of clumps, small structural differences appear between 
the clumps in Ca II H and in SJI 2796, 1330 and 1400. This multi-thermal plasma picture and 
the large density inhomogeneities at such high resolution suggests that the thermal instability 
mechanism in a strongly anisotropic magnetic held such as the corona is far more complex 
than the simple picture of a uniformly progressive cooling plasma with cool chromospheric 
cores surrounded by warmer diffuse material. On the other hand, the difference in the 
average widths between rain in emission in AIA 304 and rain in absorption in AIA 171 and 
AIA 193 suggests a relatively wide transition from transition region to coronal temperatures 


of prominences, IChiuderi fc Chiuderi Drago 


values in prominences (jVial et ah 


eri Draeo 

1991 

2012: 

Parent! 


Heinzel et al.ll200ll) . similar to reported 


20151 ). This suggests a large increase in 


complexity in the interface between chromospheric and transition region plasmas, and less 
so for the PCTR. It is also possible, although less likely, that since absorption in the AIA 171 
hlter is produced by multiple elements (but mainly H, He and He"*"), an inhomogeneous He 
distribution within the loops (more concentrated in the core of the loops than hydrogen) 
could lead to such difference in widths. How such high degree of inhomogeneity, especially 
at the smallest scales, is achieved in such loops is an important question that needs to 
be addressed in future work. Theoretically, the presence of high degree complexity at the 
smallest scales could be explained by the existence of tangential discontinuities in the held in 
which material can collapse, predicted by stability aiialysis o f prominence plasma in similar 
conditions as in the present observations flLow et al.ll2012al) . Such theoretical studies also 
predict the constant evolution of magnetic helds and huid at these small scales lead ing to 
spontaneous formation and resistive dissipation of discrete currents (ILow et al.ll2012bl) . Such 
events may explain the constant morphological a nd physical c hange s of coronal rain clumps 
that we observe along their fall, also reported in iHarra et ah! (120141) . 
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5.3. Elemental scales and global topology 

The present results strongly suggest that the bulk of the distribution is bel ow the highest 
resolution used here (~ 0.2"). This study conhrms and extends the result of IScullion et al. 


fl2014j) over various instruments spanning different temperatures. T his picture is also ob¬ 
tained in 2.5-D MHD simulations of coronal rain bv iFang et al.l (120131) . in which the number 
of rain clumps keeps increasing at seemingly ever decreasing scales. In Figs. 116] and fTTIf upper 
panel) the standard deviation of widths and the ratio between standard deviation and mean 
values for the widths decreases signihcantly from the AIA measurements down to CRISP 
measurements (except for the AIA 304 widths in the second dataset, perhaps reflecting the 
absence of large scale coronal rain downflows and no signihcant showers in this case). On 
the other hand the lengths present a more random distr ibution with peak s at a few arcsec 


or less , which is also found in numerical simulations by If^nna et al.l (1201211 and iFang et al. 
(120131) . We further hnd extremely long clumps encompassing a few tens of thousands of 
km. This large scatter in lengths may reflect the presence of several factors at play along 
the held (such as hows, thermal conduction, instabilities, held geometry etc.), supported by 


numerical results. T 


le more uniform di s tribution of widths, the agreement over a signihcant 


num 


2014 : 1 Scullion et al. 


3 er of datasets (Antolin et al.ll2012l: lAntolin fc Rouppe van der Voortll2012l: iHarra et al. 


20141) and the similarity with widths of prominence threads (ILin et al. 


20051) suggests a more fundamental mechanism behind the shaping of condensations in the 
transverse direction to the held. This behaviour strongly suggests the existence of elemen¬ 
tal strand-like structures at smaller scales. Such fundamental scales would be dehned by 
a similar thermodynamic and magnetic held evolution. A relevant question to ask is then 
whether thermal instability plays a main role in the shaping of such structures, and especially 
dehning the widths of these elemental structures, or whether the obtained widths would also 
be expected for substructure in a thermally stable loop. In Fig. [16] spatial resolution and 
spectral diagnostic (and therefore temperature) are intertwined in the results and it is dif- 
hcult to disentangle them and estimate the inhuence of each on the observed widths. The 
relation between temperature and width is shown in the lower panel of Fig. [T3 where we 
can see a roughly monotonically increasing trend. However, the fact that Ha, Ca II H and 
Mg II k form at roughly similar chromospheric temperatures but show a signihcant difference 
in widths points mostly to a major role of resolution over temperature (although less so from 
Mg II k due to the higher opacity). Also, the widths found in C II and especially in Si IV 
are similar, despite the signihcant diherence in temperature formation, suggesting that such 
small strand-like structures may extend to hotter lines at high spatial resolution. 


It is nonetheless possible that the observed strand-like structure in these loops, or part 
of it, is the product of thermal instability. For instance, in Section 14.51 we showed at the 
smallest scales we may have an inhuence of the MHD thermal mode on the shaping of strand- 
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like structure. The ripples in the clump’s neighbourhood generated by this wave imply a 
strand-like structure all around the main clump. These ripples constitute small density en¬ 
hancements (which correspond to the spatial shape of the thermal mode eigenfunctions) that 
propagate together with the clump and can enhance radiative losses. It is therefore likely 
that such p erturbations may grow into condensations themselves, leading to further strand 
generation fjSoler et al.N20121) . Besides this effect, catastrophic cooling can also introduce 
large pressure variation, especially in the strong density inhomogeneities found here with 
cool cores with electron densities on the order of 10^^ cm“^. Due to flux freezing the collapse 
of material following thermal instability generates an increase of the magnetic held. Despite 
the local character of the condensation the subsequent held enhancement is non-local due 
to the rigidity of the co ronal magnetic held (high magnetic tension due to the high external 
Alfven speed, lAnlll984j) . This process can thus generate a (slightly) magnetically diherent 
strand which would in turn shape warmer material up- and downstream. Such a scenario 
has been obtained in 2-D MH D simulations of the rmal instability in the interstellar medium, 
with similar plasma f3 values flSharma et al.l 120101) . This scenario implies a relation between 
clumps’ widths and the magnetic held at small scales, which ohers interesting coronal seis¬ 
mology applications. 

The proposed non-local enhancement of the magnetic held would increase the magnetic 
pressure downstream of the clumps, which could explain the less than free-fall speeds (and 
rather constant accelerations) always observed in coronal rain. Let us make an order of 
magnitude estimate of this possibility and leave a detailed analysis for future work. Let 
us consider the case of dataset 1 for which we have density estimates. We have an average 
velocity at the footpoint of 100 km s“^, a loop height of 42 Mm, and let us assume a negligible 
initial speed at the apex of the loop. The average ehective gravity along an elliptic loop with 
a ratio of loop height to half baseline be tween 0.5 and 2 varies between 0.132 km s“^ and 
0.216 km s“^ flAntolin fc Verwichtell20111) . Taking an average of 0.174 km s“^, the average 
velocity at the footpoint for our clumps would be around 150 km s“^ due to the effective 
gravity alone. From Bernouilli’s equation we can then estimate the increase of pressure 
due to the magnetic held downstream from the clump that would lead to a decrease of 
the hnal velocity of the clump (assuming no change in the gas pressure downstream of the 
clump). We obtain values of 0.3 — 0.5 dyn cm“^ for average downstream plasma electron 
densities of 3 — 5 x 10® cm“^ (a factor of 10 lower than the average density derived from 
the EUV absorptio n measurement in Section 14.31 and consistent with numerical simulations, 
Antolin et al.ll2010fl . These values correspond to a magnetic held increase of 2.8 — 3.6 G. On 
the other hand, the gas pressure created by the dense clump is ~ 0.38 dyn cm“^ (assuming a 
low ionisation degree), matching the obtained values and therefore supporting this possibility. 


Besides the interesting possibilities of coronal rain as local probes of the internal ther- 
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modynamic and magnetic loop structure, it can also serve as a global tracer of the magnetic 
field topology. This was hrst demonstrated in Paper 1, where the angle of fall of rain clumps 
was calculated based on the full velocity vector, which allowed to trace the lower legs of 
loops and thus to distinguish several loop families. Here we show that the rain clumps’ 
lengths may exist over a few tens of thousands of megameters, allowing fur ther tracing of 


the lo op legs. Such long rain structures have also been found recently by IScullion et ah 


(1201411 in the context of post-flare loops, forming from the apex of loops. But even if short, 
by calculating the variance over relatively long periods of time as in Fig. [TTl or by detecting 
all the clumps crossing along perpendicular cuts to the trajectory at different heights as in 
Fig. m the global shape of the loops’ magnetic held can be inferred. By doing so we detect 
a clear expansion in the loops of dataset 1 with CRISP. The width of the loop cross-section 
over which clumps are detected for loops 1 and 2 in Fig. E] decreases from 4.6" to 2.4" and 
6.8" to 3.4", respectively, from height C down to F separated over a distance of roughly 26" 
(similar to both loops). The width of the loops at the apex, detected with AIA in datasets 
1 and 2 vary around 5", suggesting that the expansions we observe with CRISP are most of 
those experienced by these loops. On the other hand, in Paper 1 little expansion is observed, 
perhaps due to the difference in the height at which rain is detected. Here we detect it down 
to chromospheric heights, just before impact, while in the previous work only detection in 
the corona was possible. The expansion seems to occur mostly uniformly between heights C 
and F, with the expansion occurring mostly above an observed kink in the trajectory of the 
clump, a little above height F. In Fig. [5l this kink is observed at y ~ 35" (for loop 2 this kink 
is at a similar height). This kink for loops 1 and 2 is visible whenever a clump is observed 
and can be followed all the way down, suggesting that it is a feature of the magnetic held 
topology above the sunspot. Similar kinks in the held above sunspots have been detected 
by measuring the change with height of the dominant period s of running penumbral waves 
from the photosphere into the chromosphere (I Jess et ahl 1201311 . 


5.4. Persistent high-speed downflows above sunspots 

The chromospheric material observed along these thermally unstable loops occurs in 
clumps and appears intermittently at coronal heights. However, this downflow has a more 
persistent and constant character at chromospheric heights, just before impact. This is 
contrary to expected for the following reasons. At low chromospheric heights projection 
ehects are stronger and visualisation of clumps relies on absorption and therefore on high 
opacity. Unless the density of clumps is signiflcant the resulting absorption would not be 
signiflcant enough for allowing detection, which then rules out a signiflcant number of clumps. 
Also, unless the Alters have sufficiently narrow spectral bands the contribution from other 
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chromospheric and photospheric emission will strongly reduce the contrast produced by the 
clumps. For these reasons, while off-limb detection of clumps can be performed in multiple 
wavelengths, on-disc detection is strongly limited. Even in the relatively high Ha opacity 
on-disc detect i on is usually strongly reduced to about 10% of the rain quantity, as shown in 
Antolin et ahl (1201211 and Paper 1. Here we show that under specihc conditions of relatively 
high Doppler velocities from the Ha clumps with respect to the background (and therefore 
reduced projection effects in the wing of Ha) detection of clumps just before impact is 
possible and the occurrence frequency of such downflows increases considerably, enough to 
give a persistent character. 

The increase in coronal rain quantity at low heights may be the combination of two 
effects. The hrst effect is the continuous cooling, achieved in the case of complete thermal 
instability down to chromospheric temperatures. As discussed previously, in this study we 
obtain evidence for continuous cooling but in a strongly inhomogeneous way in which only 
small regions of the loops cool down catastrophically. Consequently, more rain in Ha will 
be observed at lower heights. The second effect may be more important, and is strongly 
dependent on the magnetic held geometry. The increase in coronal rain quantity may be 
caused by a funnel effect due to the increase of magnetic held strength at lower heights. As 
discussed in Section 15.31 here we observe a clear expansion of about a factor of 2 in radius for 
loops 1 and 2 from heights F to C. On the other hand, only a slight change in the width of the 
clumps between both heights is found. This is shown in the upper panel of Fig. [171 where the 
average of the clumps’ widths in Ha for both loops changes from 0.305" to 0.296" from cut 
C to F (shown in Fig. [5]), with nonetheless a signihcant decrease of the standard deviation 
from 0.115" to 0.083" respectively. From conservation of mass we would then expect both an 
increase in the frequency of rain crossing the cuts at chromospheric heights and, correspond¬ 
ingly, an increase of the clumps’ lengths with decreasing height. Unfortunately, for the latter, 
we are unable to provide precise measurements of lengths at lower heights, due to the large 
projection effects and large speeds of the clumps (as shown in Fig. [T3]the height over which 
the clumps are significantly distinct from the background is limited to the low atmosphere 
and the region above the spicules). Although this tendency of elongation of clumps has been 
previously reported flAntolin et al.ll2010l. 120121: lAntolin fc Rouppe van der Voortll2012l) and 
is observed in general for dataset 2, this quantity is expected to be severely affected by other 
factors apart from the magnetic field expansion, such as flows, the effective gravity due to 
curvature, etc. This is reflected in the large scatter in lengths found for all instruments, 
shown in Fig. [161 From this funnel effect one would also expect a decrease of the falling 
speeds from the compression of plasma downstream of the blob, as predicted by numerical 
simulations (lAntolin et al.ll2010l: iMackav fc Galsgaardll200ll) . In our observations we could 
not detect a net decrease of falling speeds just before impact. However, it is possible that a 
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velocity decrease is present but cannot be confidently detected due the elongation of clumps 
and the short distance above the chromosphere over which they offer enough contrast (which 
leads to a larger uncertainty in the measurements of projected speeds). 

Persistent red sh ifts above sunspots have often been observed in transition region and 


chrom ospheric lines (iBrekke et ahl Il990l: iBruecknerl Il98ll: iDerel Il982l: iKieldseth-Moe et ah 


19881 ). often associated to bright fan-shaped structures observed i n EUV above the umbra, 
usually termed plumes f Foukal et ah 1974 : Brosius fc White 2004 : BrosiusI 2005 1. The range 
of downflow speeds associated to such structures is large, mai nly subsonic (20 — 4 0 km s~^). 


with supersonic s econd compon e nts of 100 km s ^ or more (IBrekke et ah Il990l: iBrueckner 


19811: IDerel 119821: iGurmanI 119931: iKieldseth-Moe et ahl 119881: iNicolas et ahl 1198211 . Recently, 


similar events have been detected in tran sition region lines with IRIS above the umbra of 
a sunspot, reaching speeds of 200 km s“^ flKleint et ahl 120141) . Such events, especially those 
presenting high speed components, have been mainly related to either siphon flows (resulting 
from pressure imbalance at the footpoints) or condensations flows (issued from thermal insta¬ 
bility in the coronal part of the loops). The latter scenario has been more frequently invoked 
over the years, es pecially when viewed in chromospheric lines, and since it can also gener¬ 
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At high resolution, especially in chromospheric lines, coronal rain downflows have been 
shown to have a clumpy character, thus suggesting a rather intermittent and bursty down¬ 
flow character, even if frequent, at odds with some observations of persistent and continu¬ 
ous downflows above sunspots. The present observations suggest that the intermittent and 
clumpy character of cool coronal rain (observed in chromospheric lines) can become per¬ 
sistent and continuous at chromospheric heights, thereby providing an answer to the long 
standing problem of downflows. When viewed in warmer transition region lines this char¬ 
acter may be further enhanced. Indeed, we find that in general the widths of the clumps 
increase in transition region lines, which may lead to the same effect. In AIA the darkening 
of the EUV loops is observed to have a persistent a nd continuous character. The persistency 
of the rain recently reported in iKleint et ahl (120141) also supports this scenario. 
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5.5. Role in the chromosphere-corona mass cycle 


Here we show that clumps are multi-thermal, suggesting a highly inhomogeneous density 
structure. Although a thorough statistical study of densities of coronal rain clumps is needed, 
we can make a rough estimate of the mass density flux involved in the current process for 
dataset 1. For this we assume that the mass of each conden sation is roughly p roportional to 
its length, a result supported by 3D numerical calculations flLuna et ahll2012[) . The electron 
density for a particular clump is estimated to vary between 2 x 10^° cm“^ and 2.5 x 10^^ cm“^. 
Taking an average of 1.4 x 10^^ cm“^ for all Ha clumps, an average width of 0.4", an average 
length of 3.8" and taking into account that the total number of Ha clumps is 479 for the 
two loops over a time of 95 min, we obtain an average downward mass flux for one loop of 
1.23 X 10® g s“^. For dataset 2, assuming a similar density and using the corresponding values 
found for that dataset (cf. Table [T]), we hnd a value of 5.22 x 10® g s“^. These values are on 
the same order as that estimated in Paper 1, where, the observed active region comprised 
multiple rain loops (5 or 6 loops can be distiguished, as shown in their Figs. 14 and 17) and 
a corresponding higher number of clumps (the assumed average rain density in that paper 
was however assumed 2.5 times lower than that estimat ed here). These are signihcant value s 
comparable to those found for prominence drainage (ILiu et al.l 120121: iBerger et al.ll2012l) . 
This agreement over different datasets strongly suggests an important role of coronal rain in 
the chromosphere-corona mass cycle. 


As discussed in Section [5731 the detected number of clumps in Ha may very well just be 
the tip of the iceberg, and many more unresolved clumps may exist. Let’s assume nonetheless 
that we detect most of the rain in both loops over the observed time sequence, which would 
then correspond to the amount of material in these loops involved in one limit cycle (assuming 
that most of the material in one loop is thermally unstable). We can then calculate the 
average density in these loops. For this we assume the loop to be a straight cylinder with 
an average width of 5.7" (the average between both loops of the transverse distance over 
which clumps are detected at heights C) and a loop length of 130 Mm, as estimated in 
Section 13.21 (assuming a similar length for both loops). We obtain an average electron 
density of 6.2 x 10® cm“®. Using corresponding values for the loop in dataset 2 a value of 
7.3 X 10® cm“® is obtained. These values are at least 5 times below the expected density of 
active region loops in a thermal non-equilibrium state, obtained from numerical simulations. 
This suggests three possibilities. In the hrst and second the loop does not become fully 
thermally unstable with only parts of it condensing and falling as coronal rain. In the hrst 
scenario, most of the loop’s mass remains at coronal temperatures and undergoes a different 
thermodynamical evolution. In the second, the entire loop cools but only a small part of it 
reaches chromospheric temperatures while most of it stays at transition region temperatures. 
The third scenario consists of a fully thermally unstable loop, with most of the mass reaching 
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chromospheric temperatures but iuvolviug raiu clumps smaller thau those preseutly resolved. 
The first sceuario would imply that ueighbouriug strauds do uot geuerally cool dowu at the 
same time, somethiug at odds with previous (Paper 1) aud the preseut results (especially for 
dataset 1, while for dataset 2 it is possible that a siguihcaut amouut of material at coroual 
temperatures is still preseut). lu the secoud sceuario, we would expect the bulk of the raiu 
distributiou to cousist of less deuse aud warmer clumps beiug larger aud louger (for iustauce, 
3 times larger aud louger clumps aud twice the curreutly detected uumber would suffice). 
While this roughly matches the AIA statistics for dataset 1, it is however at odds with the 
IRIS results of dataset 2, which show a very similar morphology aud distributiou of clumps 
iu chromospheric aud trausitiou regiou temperatures, aud especially far fewer uumber of 
clumps at high temperatures. Heuce, the most likely possibility is the third, aud iuvolves a 
far larger uumber of uuresolved chromospheric clumps. This result suggests that the maiu 
coutributor to dowuflows (iu terms of mass) iu thermally uustable loops is coroual raiu. This 
sceuario should be tested iu 3D uumerical simulatious. 


6. Conclusions 

In this paper two datasets that combine multi-wavelength observations with SST /CRISP, 
Hinode/SOT, IRIS/SJl and SDO/AIA, spanning chromospheric, transition region and coro¬ 
nal temperatures, are used to reveal the multi-strand and multi-thermal aspects of coronal 
rain at unprecedented detail. We provide important characteristics of the thermally unstable 
plasma ending up as coronal rain, and more generally of the thermal instability mechanism 
in the corona. We show how this mechanism can play a siguihcaut role into the shaping 
of coronal structures. Furthermore, we make hrst steps towards estimating the fraction of 
coronal plasma in a normal active region undergoing full catastrophic cooling. The results 
can be listed as follows: 

• While the standard model of progressive cooling is observed, in which light curves peak in 
order of decreasing temperature, we also hud a case in which coronal emission is present 
and is even enhanced co-temporally and co-spatially with the chromospheric and transition 
region emission from coronal rain. Loops in a thermal non-equilibrium state may therefore 
not forcibly show gradual systematic cooling from coronal to chromospheric temperatures, 
i.e. a sequential EUV intensity variation from hot to cool temperatures (where the hot 
emission is signihcantly reduced) is not a necessary observational condition of thermally 
unstable loops. 

• EUV intensity darkening is found strongly correlated in space to coronal rain appearance. 
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Furthermore, the rain can produce EUV intensity variations on short timescales of a few 
minutes that can appear quasi-periodic over time. 

• Progressive cooling of coronal rain is observed, leading to a height dependence of the 
emission. Catastrophic cooling is observed to follow two stages: a rapid cooling throughout 
transition region temperatures and a slower cooling down to chromospheric temperatures, 
which may reflect the transition to optically thick plasma states. 

• Coronal rain is a highly multi-thermal phenomenon encompassing transition region and 
chromospheric temperatures. We hnd a high degree of co-spatiality in the multi-wavelength 
emission implying a strong degree of density inhomogeneity within thermally unstable 
loops. A thin transition from chromospheric to transition region temperatures must exist 
on spatial scales lower than 0.33". 

• We estimate coronal rain core electron densities to vary between 2 x 10^° cm“^ and 2.5 x 
10^^ cm“^. 


• The distribution of coronal rain widths is found to be rather independent of tempera¬ 
ture with a clear increase of clump number (and a decrease of standard deviation) at 
small scales for each wavelength. The number of clumps is, however, found to increase 
sharply across the temperatures towards the chromospheric range, reinforcing the tip-of- 
the-iceberg scenario proposed in Paper 1. The present results strongly suggest that the 
bulk of the distribution is below the highest resolution used here (~ 0.2"). 


The coronal rain clumps appear organised in strands. Such structure is not limited to chro¬ 
mospheric temperatures but extends at least to transition temperatures as well, suggesting 
an important role of thermal instability in the shaping of fundamental loop substructure. 
At the smallest detected scales we hnd structu re reminiscent of the MH D thermal mode 
flFieldlIl965l : Ivan der Linden fc GoossensIllOOlalibl: iMurawski et al.ll201ll) . 


• Through order of magnitude estimates we show that the strong density inhomogeneity that 
we hnd entails a local increase of gas pressure accompanied by a non-local increase of mag¬ 
netic pressure due to hux freezing all along the clump’s trajectory (up- and downstream), 
thereby shaping hotter material into strands. We further suggest that this mechanism 
may explain the commonly observed lower than free-fall speeds of coronal rain. 


We show how coronal rain can serve as a global tracer of the magnetic held topology. We 
detect the presence of a clear kink in the held in th e upper chromosp here above sunspots, 
reminiscent of the sharp change in angle reported in I Jess et ahl ( 20131) in a similar location. 
The magnetic held is found to expand roughly uniformly above this kink up to coronal 
heights over at least a factor of 2 in radius. 
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• The intermittent and clnmpy appearance of coronal rain at coronal heights becomes more 
continnous and persistent at chromospheric heights just before impact. This is likely due 
to progressive cooling of the rain and, especially, to a funnel effect from the expansion of 
the magnetic held at low heights. 

• The persistent (and burst-less) character of coronal rain at low heights is in agreement 
with the common detection of moderate and high speed red shifts above sunspots in 
chromospheric and transition region lines, and may thus offer an explanation to this long 
standing problem. 

• In terms of hlling factor (and therefore detectability, especially with poor resolution in¬ 
struments), coronal rain is mainly a transition region phenomenon. However, in terms of 
mass, the chromospheric component of the cooling plasma plays a major role in the mass 
cycle of the loop. We hnd average downward mass huxes per loop of 1 — 5 x 10® g s 

i n accordance with previous estima tes in Paper 1, and prominence mass drainage values 


(Liu et al. 

2012; 

Berber et ah 

2012) 
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Fig. 1.— Composite image from SDO/AIA combining three AIA filters, 193 {red), 171 
(green) and 304 (blue), showing dataset 1; AR 11903 at the East limb on the 26th of June 
2010. The dotted square denotes the Eeld-of-view of Fig. [3l The dotted loop corresponds to 
loop 2 studied in Section 13.21 The dashed square denotes the held-of-view of the SST. An 
animation of this figure is available in the online material as Movie 1. 
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Fig. 2.— Composite image combining AIA 171 (blue), IRIS/SJI 1400 (green) and Hin- 
ode/SOT Ca II H (red), showing dataset 2: AR 11903 at the West limb on the 29th Novem¬ 
ber 2013. The dashed square denotes the held-of-view covered by Hinode/SOT. The loop 
with coronal rain towards the centre right of the square corresponds to the loop studied in 
this dataset. An animation of this hgure is available in the online material as Movie 2. 
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Fig. 3.— Evolution of a loop hosting coronal rain in three AIA diagnostics and CRISP Ha for 
dataset 1. The top two panels show the loop (dehned as loop 2 in the text) at two instances 
in time (indicated by vertical dashed lines in the lower panel) in AIA 193 (red), 171 {green), 
and 304 {blue), with the loop path traced by the curved dashed line in the upper two panels. 
The long-dashed box outlines the extent of the SST £eld-of-view, the lower left corner of 
which coincides with that of the panels. The lower panel shows the counts normalised to 
their average as function of time for AIA 193 {long-dashed red), 171 {solid green) and 304 
{dashed blue), where these intensities have been obtained over the box around the loop top 
in the upper panels. For CRISP Ha {solid black) we plot a histogram of intensities (scaled 
by 0.5 and shifted by -1-0.8 to fit nicely in the panel) that stem from measurements across 
the short white line perpendicular to the loop path located at about (155", 100") (which 
corresponds to cut C2 in Fig. [5]), and are weighted by the respective clumps’ widths squared 
(see Fig. [7] for further details). 
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Fig. 4.— Multi-panel figure showing a snapshot of a clump falling towards the sunspot in 
dataset 1. From left to right, a Doppler image from SST /CRISP in Ha ±0.76 A, SDO/AIA 
304, 171, 193, 211 and 335. The dotted lines follow the centre of the condensation, observed 
in emission above the spicular layer (above y ~ 47") and in absorption below this layer. 
Notice the kink in the trajectory of the clump at ~ 35". 
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Fig. 5.— Composite image combining a Doppler image from SST /CRISP in Ha ± 0.76 A 
(in red), SDO/AIA 171 (in green) and SDO/A\A 304 (in blue) for dataset 1. The image 
shows the same snapshot of the falling clump as in Fig. 01 whose trajectory is indicated 
by the dotted curve. Four transverse cuts to the main axis of two rainy loops are shown, 
located at two different heights: Cl and FI for loop 1 (C2 and F2 for loop 2) at coronal and 
chromospheric heights, respectively. The observed clump belongs to loop 1. Movies 3 — 6 
are zooms into this figure, showing the detected clumps crossing cuts Cl, C2, FI and F2 
respectively (where the FWHM of the fit to the clumps’ widths is plotted over in red). 
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Fig. 6.— Histograms of temperature (left) and Doppler velocity (right) for the Ha clumps 
detected with SST/CRISP in dataset 1, following the method described in Section [231 
The black histogram corresponds to events whose measurements have a standard deviation 
above 10% of the mean. Positive and negative Doppler velocities correspond, respectively, 
to redshifts and blueshifts. 












46 


AIA 304 - AIA 193 - Ha 

AIA 131 - AIA 211 

AIA 171 - AIA 335 



0 20 40 60 80 

Time [min] 


Fig. 7.— Time evolution of the intensity across cuts Cl and C2 (corresponding to loops 1 and 
2 in dataset 1, respectively) for multiple AIA filters (indicated at the top with their respective 
colour) and the condensations detected with CRISP in Ha (in black). The histogram of Ha 
intensities is constructed by taking a bin size of a 100 s, and for each time bin the intensity 
of all contained clumps are summed, each multiplied by the clump’s width squared. The 
AIA intensities at each time step are averages over a region of 1.8" in width in the transverse 
direction to the cuts and same length as the cut, except at times where condensations are 
detected in Ha. For these, the average is not calculated over the entire cut but limited to the 
positions along the cut where the clump is detected (dehned by the width and the location 
of the clump). 
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Fig. 8.— Multi-panel showing a snapshot of the coronal rain event in dataset 2. From left to 
right we have ifmode/SOT in Ca II H, IRIS/SJI 2796, 1330, 1400 and SDO/AIA 304. Each 
snapshot corresponds to the closest snapshot in time to that taken by SOT, and for each 
the average background is subtracted. The time difference of each snapshot with respect to 
that of SOT is indicated in the upper part of each panel. An animation of this hgure can be 
found in the online material as Movie 7. 
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Fig. 9.— Light curves {top panel) and cumulative light curves {bottom panel) of the coronal 
rain pixels in the loop highlighted by the rain in Fig. [8] of dataset 2. In both panels the 
Hinode, IRIS and SDO diagnostics are differentiated by colour: Ca II H {purple), C II 1330 
slit-jaw {blue), Si IV 1400 slit-jaw {cyan). Mg II k 2796 slit-jaw {green), AIA 304 {yellow), 
and AIA 171 {red). The SOT and AIA light curves have been smoothed over 7 and 3 time 
steps, respectively, close to the cadence of the IRIS slit-jaw images. 
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Fig. 10.— Composite image of the rain event in dataset 2 combining i/mode/SOT in Ca II H 
(in red), /ii/^/SJI 1400 (in green) and IRIS/SJI 1330 (in blue). An animation of this hgure 
can be found in the online material as Movie 8. 
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Fig. 11.— Substructure for a rain shower in the rain event of dataset 2. The shower 
corresponds to the head of the rain and is shown at three different times separated by 
roughly 38 s from each other (shown in the rows) and for different wavelengths (shown in 
the columns). From left to right we have hfmode/SOT in Ca II H, IRIS/S31 2796, 1330, 
1400 and SDO/AIA 304. The SJI and AIA snapshots correspond to those closest in time to 
the SOT snapshot in the hrst column of each row. The time difference with respect to the 
hrst snapshot (top left panel) is indicated in the bottom of each panel. 
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Fig. 12.— The L function (in logarithmic scale), as dehned in iLandi fc Realel (120131 ). with 
respect to temperature, for the various AIA hlters calculated over the trajectory of a rain 
condensation producing EUV darkening. 
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Fig. 13.— Histograms of rain condensations detected in Ha with SST /CRISP for dataset 1 
along cuts Cl, FI and C2, F2 for loops 1 and 2, respectively. The condensations detected 
at coronal heights (C cuts) constitute the deep blue histogram, while those detected at 
chromospheric heights (F cuts) constitute the light blue histogram. The histograms are 
constructed by taking a bin size of a 100 s, and for each time bin the intensity of all contained 
clumps are summed, each multiplied by the clump’s width squared. 
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Fig. 14.— Variance of a Doppler image in Ha ± 0.6 A with ^FT/CRISP over a time of 
22 min for loop 1. Varions strands are visible from 50"to 55" and 28" to 35" in the vertical 


axis. 
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Fig. 15.— A rain condensation captured in Ha by SST/CKLSP, here shown as a Doppler 
image {left panel) at ±1.1 A. The intensity of the image is saturated in order to better 
visualise the small dim structure parallel to the main condensation. Transverse cuts to the 
direction of propagation of the condensation are shown towards the head (location A) and 
tail (location B). The intensity along these cuts is shown in the right panels at both A and 
B locations as stack curves separated by an ad-hoc value. The colour of the intensity curve 
corresponds to the cut of the same colour in the left panel (for location A or B). 























- 55 - 


>^ 

u 

c 

CD 

D 

CT 

(L) 


<U 

> 

o 

0 ) 

cr 


0.5 

0.4 

0.3 

0.2 

0.1 

0.0 



- CRISP Ha - 

- SOT Coll H - 

SJI/2796 

- SJI/1330 

- SJI/1400 “ 


0.30 

0.25 

0.20 

0.15 

0.10 

0.05 

0.00 



0.0 0.5 1.0 1.5 


0.0 0.5 1.0 1.5 2.0 


Width [orcsec] 


Length [orcsec] 




0.02 

0.00 


0 


10 


15 20 25 30 


Width [orcsec] 


Length [orcsec] 


Fig. 16.— Histograms of calculated widths (left) and lengths {right) for the detected coronal 
rain condensations in the various wavelengths used in the present study for datasets 1 and 
2. The corresponding colour for each wavelength is shown in the legend. The bin size in 
each wavelength is set to the spatial resolution of each instrument. For the high resolution 
instruments {upper panels) only the small scale part of the spatial range is shown. In the 
lower panels only the AIA results are shown. 
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Fig. 17.— Diagrams showing relationship between various quantities. As quantities we have 
the height of the measurements, estimated from the footpoint location of the loop in the 
chromosphere, the width of the condensations and the temperature of the condensations. 
For each quantity Q its average value (Q) over the entire set of measurements is plotted, 
together with error bars corresponding to the standard deviation a {{Q) ± a). Each colour 
corresponds to a different wavelength, indicated in the top panel. 





































